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ABSTRACT 



We present the analysis of the galaxy stellar mass function in different environments at intermediate redshift (0.3 < z < 0.8) for two 
mass-limited galaxy samples. We use the IMACS Cluster Building Survey (ICBS), at masses M, > 10'""^ Mq, to study cluster, group, 
and field galaxies at z = 0.3 - 0.45, and the ESO Distant Cluster Survey (EDisCS), at masses M, > lO'^-^Mg, to investigate cluster 
and group galaxies at z = 0.4 - 0.8. Therefore, in our analysis we include galaxies that are slightly less massive than the Milky Way. 
Having excluded the brightest cluster galaxies, we show that the shape of the mass distribution does not seem to depend on global 
environment. Our two main results are: (1) Galaxies in the virialized regions of clusters, in groups, and in the field follow a similar 
mass distribution. (2) Comparing both ICBS and EDisCS mass functions to mass functions in the local Universe, we find evolution 
from z ~ 0.4 - 0.6 to z ~ 0.07. The population of low-mass galaxies has proportionally grown with time with respect to that of massive 
galaxies. This evolution is independent of environment ~ the same for clusters and the field. Furthermore, considering only clusters, 
we find that no differences can be detected neither within the virialized regions, nor when we compare galaxies within and outside 
the virial radius. Subdividing galaxies in terms of color, in clusters, groups, and field red and blue galaxies are regulated by different 
mass functions, but comparing separately the blue and red mass functions in different environments, no differences are detected in 
their shape. 

Key words, galaxies: clusters: general — galaxies: evolution — galaxies: formation — galaxies:luminosity function, mass function 



1. Introduction 

In standard A cold dark matter (ACDM) cosmological mod- 
els, cold dark matter haloes form from the gravitational col- 
lapse of dark matter around peaks in the initial density field. 
Haloes assemble hierarchically, such that smaller haloes merge 
to form larger and more massive haloes in d ense environments 
(iMo & Whitelll996l: ISheth & Tormenll2002h . According to the 
current paradigm of galaxy formation, galaxies form within 
haloes, owing to the cooling of hot gas. Haloes and galaxies 
evolve simultaneously, and the evolution of a galaxy is driven by 
the evolution of its host halo. If the halo is accreted by a larger 
halo, the galaxy will be affected by it as well: for example, the 
galaxy's diffuse hot gas reservoir may be stripped, removing its 
fuel for future star formation (e.g. [La rson, Tinslev, & Caldwell 
1980; Balogh. Navarro. & Morriil200Q[ iWeinmann et alJl2006t 
van den Bosch et al. 2008). The evolution can also be governed 
by the interplay between smooth and clumpy cold streams, 
disk instability, and bulge formation. Intense, relatively smooth 
streams maintain an unstable dense gas-rich disk. Instability 
with high turbulence and giant clumps is self-regul ated by 
gravitational interactions within the disk (see, e.g., Dekel et alJ 



l2009l and references therein). Galaxies may also experience ma- 
jor mergers, which transform late-type gala xies into early-type 
galaxies with a central bulge component (e.g. lDriver et alj|2()06l : 
iDrorv & Fisheil l2007h . Mergers drive gas towards the centre, 
where it can trigger a burst of star formation and fuel the central 
black hole, the feedback from which can he at the remaining gas 



and e ventually quench s tar formation (e.g. [Mihos & H ernauist 



IWild et alJ l2507l: iPasquali et alJ l2()08l: Ischawinski et al 
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Several studies have shown that there are also other exter- 
nal stresses that act mainly on galaxies in dense environments 
and, in general, do not allow the maintenance of spiral structu re. 
For example, ram pressure dGunn & Gottll 19721 : iBek ki 2009) is 
a drag force that is capable of stripping the galaxy of much of 
its interstellar gas and preventing the formation of new stars . 
Galaxy harassment (iMoore et al.ll996tlBoselli & Gavazzil2006l) 
is a mechanism that strips a galaxy of part of its mass and drives 
a morphological transformation as a consequence of frequent 
high speed encounters. Harassment has the potential to change 
any internal property of a galaxy within a cluster, including the 
gas distribution and content, the orbital distribution of stars, and 
the overall shape. Finally, cluster tidal forces (Bvrd & Valton^ 
[l990) can act with different efficiency depending on environ- 
ment, such that it could be possible that field galaxies infalling 
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into larger structures can be transformed for example from gas- 
rich spirals to gas-poor lenticular galaxies. 

Hence, in their evolution, galaxies are expected to be 
strongly influenced by the environment in which they reside. 
Several works have asserted that in their evolution galaxies 
are also and mostly influenced by their stellar mass. For ex- 
ample, Kauffmann et al.l ( l2003h found that color, specific star 
formation rate, and inter nal structure are stron gly correlated 
with galaxy stellar mass. iPasquah et al.l (l2009l) demonstrated 
that the star formation and AGN activity of galaxies have a 
much stronge r dep endence on stellar mass than on halo mass. 
iThomas et alJ (1201 0) argued that the formation of early-type 
galaxies is environment-independent and driven only by self- 
regulation processes and intrinsic galaxy properties such as 
mass. 

Distinguishing the separate contributions of environment 
processes with those driven by an intrinsic property is clearly 
critical to understanding galaxy evolution. In the nature versus 
nurture scenario, mass represents the primary "intrinsic prop- 
erty" closely related to primordial conditions, while the environ- 
ment represents all the several possible external processes that 
can influence galaxies in their evolution. 

In the literature, there are many works in which global en- 
vironme nt and mass are segregaited and analysed separately 
(see e.g. van den Bosch et al. 2008; Guo et al. 2009'.' lovino et alJ 
12010; Mercurio et al. 2010; Peng et al. 2010), to study die vari- 
ations in galaxy properties as a function of either mass or envi- 
ronment, by fixing one of these and studying how the properties 
vary as the other changes. However, studies of galaxy properties 
as a function of mass normally have not considered the possibil- 
ity that the galaxy mass distribution itself may vary with envi- 
ronment. 

In this work, we investigate whether and how the effect of 
mass and environment are related, and whether environment can 
influence masses, in particular the galaxy stellar mass distribu- 
tion. While some estimate of galaxy stellar mass can be quite 
easily obtained (even though uncertainties are large), there are 
different ways to describe the environment. For masses, it has 
been shown that determinations resulti ng from different met hods 
agree well within the errors (see, e.g. lBeU & de Jongll200ll who 
give a relation bet ween a stellar mass-lu minosity ratio and the 
galaxy colour, and lBolzonella et ani2010l who use a SED fitting 
technique using the code Hyperzmass, a modified version of the 
photo metric redshift code Hyperz ( Bolzonella. Miralles, & Pelldl 
I2000h ). This is true, although different choices of IMF, stel- 
lar population synthesis model, SFR history, metalhcity and 
adopted extinction curve (see, e.g., iMarchesini et al.l l2009t 
iMuzzin et al] l2009'; C onrov et al.1 l20l0') can be expected to sys- 
tematically affect the mass estimates. For the environment, it is 
possible to refer to e ither the globa l or loc al environment. As 
discussed in detail in ' Muldrew et al.l (1201 l l). there is no univer- 
sal environment measure and the most suitable method depends 
on the scale being probed. In the case of the global environment, 
galaxies are commonly subdivided into e.g. superclusters, clus- 
ters, groups, field galaxies, and voids, according to the host halo 
mass, while in the case of the local environment, it is described 
through the estimates of the local density, which can be calcu- 
lated following several definitions. 

Several works have focused on the galaxy mass distribu- 
tion and its evolution in one global environment - the field, 
but very little is known about the mass function in clus- 
ters of galaxies. Studies focused only on the field, such as 
[Drorvetal. (2005); Gwvn & Hartwick (2005); Fontanaetal. 
(>2006.) ; .Bundv et al.. (.2006.) : Pozzetti et al.i (.2007,) : .Drorv et ah 



(l2009h . and iBaldrv et al] (l2012h . presented the mass func- 
tion of all galaxies and their results are in good agreement. 
Fontanaetal. (2004, 2006); Bundy et al. (2006); Borch et al| 
(2006) and iPozzetti et all (120101) demonstrated that for galax- 
ies with M, > IO'^Mq, overall, the evolution of the total mass 
function from z - 1 to z = is relatively modest, which im- 
plies that the evolution of objects with mass close to the local 
characteristic mass is essentially complete by z~l. On the other 
hand, they found that less massive galaxies evolve more than 
massive ones, displaying a rapid rise beyond z~l (for the mass 
function at z > 1 see Eisner et al. 2008; Kajisawa et al. 20(39 



Marchesini et al. 2( )09l;ICaputi et al.i201 ItlGonzalez et al.li201 1 
Mortlock et al...201 i1k 



Some works have also analyzed galaxies of different mor- 
phological types separately, since it is known that galaxies of 
different types and with different star formation histories con- 
tribute in different ways to the mass function, shaping, for ex- 
ample, either the massive tail or the low mass end of the total 
mass distribution. There are several ways to subdivide galaxies 
into at least two populations (i.e. early and late-types), usually 
either according to (1) their star formation histories (being pas- 
sive or star-forming, for example using a rest-frame color, hence 
separating blue and red galaxies, their SEDs, their spectroscopic 
features), or (2) their structure (structural parameters or mor- 
phologies). [Baloghetal] (12001') analyzed the environmental de- 
pendence of the luminosity function and the associated stellar 
mass function of passive and star forming galaxies in the Two 
Micron All Sky Survey, and found that, in the field, active galax- 
ies follow a much steeper high mass end mass function than pas- 
sive galaxies, while in clusters both active and passive galax- 
ies have a steep high mass end. Sub dividi ng field g alaxies de- 
pen ding on flieir colors.lBaldrv et al.i(l2004 : lBaidrve t al. (200^ 
and EaWrveTaD (12001 in the local Universe, and feorch et all 
(200l) and iBolzonellagtaD (1201 Oh at intermediate redshifts, 
identified a bimodal shape in the mass function in the field 
with an upturn related to the two different populations: early- 
type galaxies dominate the high masses, while late-type galax- 
ies mostly contribute to the intermediate/low-mass part of the 
mass function at all redshifts. Moreover, they have shown that 
the mass functions of early-type and late-type galaxies evolve 
differently with redshift. 



As far as clusters are concerned, in IVulcanietal](l2011h for 
the first time we studied the mass function in clusters, find- 
ing a quite strong evolution with redshift. Clusters in the lo- 
cal Universe are proportionally more populated by low mass 
galaxies than clusters at high z. We concluded that (1) the mass 
growth caused by star formation plays a crucial role in driving 
the evolution; (2) it must be accompanied by infall of galaxies 
onto clusters; and (3) the mass distribution of infalling galax- 
ies might be different from that of cluster galaxies. To this aim, 
we compared our results for clusters with the field mass func- 
tions found in the literature, to see whether galaxies in differ- 
ent global environment are characterized by different mass func- 
tions. In our preliminary analysis, we found that at high masses 
(log M^/Mq > 11), the mass functions of field and cluster galax- 
ies at high-z have rather similar shapes, while the situation at 
intermediate-to- low masses is ambiguous. Indeed we found that 
different field studies give quite different results at these masses. 
If we followed Ilbertet al. (2010), we could suggest that field 
galaxies have a steeper mass function than cluster galaxies, indi- 
cating the presence of a significant environmental mass segrega- 



2 



Benedetta Vulcani et al.: No dependance of the mass function on global environment 



tion. In contrast, the results of Bundy (200513 suggest that there 
are no large differences between the mass distribution of galaxies 
in the different environments at high-z. Unfortunately, based on 
these results, it remains unclear whether field and cluster galax- 
ies have similar or different mass distrib utions. The preliminary 
results presented in IVulcani et alj (1201 Ih were performed using 
inhomogeneous data and slightly different redshift ranges, so it 
can not be used to draw definite conclusions. 

On the theoretical side, Mo ster et al.l ([2010) found a corre- 
lation between the stellar mass of the central galaxy and the 
mass of the dark matter halo. Using N-body simulations, they 
found that the clustering properties of galaxies are predomi- 
nantly driven by the clustering of the halos and subhaloes in 
which they reside, and provided a model to predict clustering 
as a function of stellar mass at any redshift. This result could 
also suggest that also the total (central+ satellites) mass function 
may depend on environment. However, the correlation between 
the total galaxy stellar mass function and the mass of the par- 
ent halo has not yet been studied. It would be very interesting to 
understand whether simulations predict a mass segregation with 
the environment, considering the initial and evolved halo mass 
and how they predict the evolution with redshift as a function 
of the environment. This would allow us to understand the role 
of the mass halo in influencing the evolution of galaxy masses. 
This analysis is deferred to a forthcoming paper (Vulcani et al., 
in preparation). 

2. Aims of this work 

So far, observational studies have shown that in different en- 
vironments early and late-type galaxies follow different mass 
distributions and they are also present in different proportions. 
As mentioned above, the theoretical expectations suggest that 
galaxy mass of the central galaxy depends on environment. This 
may have raised to the expectations that also the total mass func- 
tion is different in different environments. The main goal of our 
work is to test this proposition and study the stellar mass distri- 
bution as a function of the halo mass: indeed, to compare mass 
functions of galaxies in clusters and field means also to compare 
mass functions of galaxies hosted in haloes of different masses. 

We want to investigate whether the mass function is "uni- 
versal" and, if this is the case, how this comes about. The main 
questions of this paper can be summarized in this way: do obser- 
vations suggest that the mass function at intermediate redshifts is 
driven by the halo mass? Does the mass function of red and blue 
galaxies separately depend on halo mass? Does the evolution of 
the mass functions depend on global environment, therefore is it 
expected to depend on halo mass? Or, as an alternative, is it pos- 
sible that the galaxy mass distribution is unaffected by where it 
has formed and which halo it comes from? Is there some mech- 
anism that gives the same imprinting to all galaxies, regardless 
where they are? 

The analysis p resented in thi s work is complementary to the 
work presented in I Vulcani et al.l ( 12012 '). There, we analyze the 
role of local density in shaping the mass function, using a near- 
est neighbour-based measure, that is largely in dependent of dark 
matter halo mass, as shown bv lMuldrew et alj ( |201 1 ). In that pa- 
per, we address the following questions: does the mass function 
depend on local density at low- and intermediate-z? In both the 
field and clusters? How does the mass function change with local 
density? 

' From private c om munication, these d ata are the combination of 
iBundv et alJ (l2005l) and lBundv etll\ ( l2006h . 



The two papers therefore address different points. The two 
ways to define the "environment" are not equivalent, and in fact 
are giving different information (halo masses versus local phe- 
nomena). Sometimes the differences between local and global 
environment are subtle and confusing, so it is possible that part 
of the results presented in this paper will surprise the reader. 
They are not in line with what it is generally expected. We will 
show in this paper that, even if the mass function d oes depend 
on local density (as we show in lVulcani et al.ll2"012l) . the differ- 
ences in local density distributions in clusters compared to the 
field are insufficient to induce a difference in the mass functions 
in these global environments. Thus, the investigation of the de- 
pendence of mass functions on global environment that we make 
in this paper is an independent test of whether the global envi- 
ronment al one is able to produce a difference in the galaxy mass 
function. In lVulcani et alJ (1201 2h we also put together the results 
of that and this works, contrasting the role of global and local 
environments in shaping the mass functions. 

The study of the mass functions has been developed only in 
the last years, while much more effort has been spent to charac- 
terize the luminosity function. In the literature, there are several 
works that carefully analyzed luminosity functions, in terms of 
both the evolution with redshift and the dependence on the envi- 
ronment. It might be reasonable to expect that the mass function 
is simply a mirror of the luminosity function, being the mass 
strictly linked to the luminosity of the galaxy, hence it is com- 
monly thought that results on the luminosity function can be 
used to infer also the mass function. On the contrary, there is 
not a linear correlation between mass and luminosity: galaxies 
do not all have the same colors, hence one single mass-to-light 
ratio, as may be true for passively evolving galaxies. The con- 
sequence is that the luminosity function does not provide direct 
information about the mass function (see Appendix A) and they 
can also not give the same results. 

Moreover, it should be kept in mind that the samples used 
for studying luminosity and mass functions should be assembled 
following different criteria. Traditionally, luminosity functions 
are studied in magnitude limited samples, where a cut in lumi- 
nosity is performed. Instead, the best choice to study mass func- 
tion is to adopt mass limited samples, that include all galaxies 
more massive than a limit, regardless of t heir color or morpho - 
logical type. As discuss in Appendix B of IVulcani et alj (1201 lb . 
it is important to note that the choice of a magnitude limit im- 
plies a natural mass limit below which the sample is incomplete. 
Hence the mass distribution derived from a magnitude limited 
sample is meaningless, because affected by incompleteness, be- 
low the the limit corresponding to the mass of a galaxy with the 
reddest color and the faintest magnitude in the sample. 

As a consequence, the characterization of the mass distribu- 
tion in a mass limited sample is very important to study the role 
of the mass in driving several galaxy properties and to under- 
stand how much the environment can influence the mass distri- 
bution at different cosmic epochs. 

Our goal for this paper, then, is to compare the galaxy stel- 
lar mass distribution in cluster regions, cluster infalling regions, 
groups, and the field using homogeneous data, in order to es- 
tablish whether and by how much the total galaxy stellar mass 
function depends on global environment at a fixed redshift. To 
understand the evolution of the mass function in clusters, it is 
especially useful to consider the mass function in the cluster out- 
skirts, where galaxies will have time to become part of clusters 
before z = 0. In this work, we also perform a cut in color, to sep- 
arate star-forming from passive galaxies, to analyze for the first 
time the mass function of blue and red galaxies also in clusters. 
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The paper is organized as follows: in ^we present the sur- 
veys we use for the analysis — the IMACS Cluster Building 
Survey (ICBS) and the ESO Distant Cluster Survey (EDisCS) 
— and depict their main characteristics. In ^we define the en- 
vironments we analyze, in^we show our results: in ^S.ll we an- 
alyze the galaxy stellar mass function as a function of the global 
environment, in ^5.2l we compare our findings with some results 
from the literature to depict the evolution of the mass functions 
in different environments. Then, we analyze the galaxy stellar 
mass function in clusters ( 35.31 ), and as a function of the colour 
( 35.41) . In ^we discuss our results, explaining the implications 
of our finding in the evolution of mass functions ( 36.11 ), and in 
the dependence of the mass distribution on galaxy properties 
( 36.21 ). We also contrast the different role of the global and local 
environments in shaping the mass function ( 36.3l l. Finally, in ^ 
we summarize our results. 

Throughout this paper, we assume - 70kms"' Mpc"', 
Qm = 0.30, and - 0.70. The adopted initial mass function is 
that o f lKrouDal(l2001h in the mass range 0.1-100 Mq. All magni- 
tudes used in this paper are Vega magnitudes. 



3. DATA SET 

In this paper, we take advantage of two different surveys to per- 
form an analysis of the mass function. We use the ICBS data 
to characterize galaxies at intermediate redshifts (0.3 < z < 
0.45) in different environments. These data are complemented by 
EDisCS data to study a large sample of galaxies at 0.4 < z < 0.8. 

The ICBS provides homogeneous spectroscopic data of 
galaxies in several environments. In this way, the redshift mea- 
surements are very accurate, being derived from the spectro- 
scopic analysis, and the membership to the different environ- 
ments is also well established. 

EDisCS contains a much larger sample of cluster and group 
galaxies, although spectroscopic redshifts are available for only 
a subset of them. Photometric redshifts are theref ore used, even 
though they are less reliable. In Appendix A of I Vulcani et alj 
(1201 ih . we illustrated that the galaxy mass function determined 
from photo-z's and photo-z membership agrees with the mass 
function determined using only spectroscopic members and 
spectroscopic completeness weights, in the mass range across 
which they overlap. 

3.1. ICBS 

The IMACS Cluster Building Survey (ICBS) (Oemler et al. 
2012a, submitted) is focused on the study of galaxy evolution 
and infall onto clusters from a clustercentric radius R ~5 Mpc to 
the cluster inner cores. Data have been acquired using the wide 
field of the Inamori-Magellan Areal Camera and Spectrograph 
(IMACS) on Magellan-Baade for four fields containing clusters 
at intermediate redshift. 

The ICBS sought to define a homogeneous sample of clus- 
ters by selecting the most massive cluster per comoving volume 
at any redshift. Clusters were selected using the Red-Sequence 
Cluster Survey method Gladders & Yee (2000), either from the 
RCS itself, or from the Sloan Digital Sky Survey in regions 
of the sky not covered by the RCS. As described in Oemler 
et al. (2012a, submitted), direct imaging in the griz bands was 
obtained for two of the fields with the f/2 camera of IMACS. 
Imaging in the B VRJ bands was obtained for the other two fields 
using the Wide Field CCD camera on the du Pont Telescope. In 
addition, very deep r-band photometry, complete to r = 25.0 was 



obtained for all fields with IMACS. Spectroscopic targets were 
selected from r-band photometry down to a limiting magnitude 
r = 22.5 for all fields. 

The IMACS f/2 spectra have an observed-frame resolution of 
10 A full width at half-maximum with a typical S /N ~ 20 - 30 
in the continuum per resolution element. In each 28' diameter 
IMACS field, spectra for 65% of the galaxies that are brighter 
than r ~22.5 were taken on the 6.5m Baade Telescope at Las 
Campanas. Of those observed, only about 20% failed to yield 
redshifts, or turned out to be stars. 

Details of the data and its analysis are presented in Oemler 
et al. (2012a, submitted ) and Oemler et al. (2012b, submitted). 

The data discussed in this paper come from four fields that 
contain rich galaxy clusters at z = 0.33, 0.38, 0.42, and 0.43, 
as well as other structures at different redshifts. In this paper, 
we decided to restrict our analysis to ICBS galaxies in the red- 
shift range 0.3 < z < 0.45, in all the environments treated. This 
was done to focus on a rather limited redshift range in order 
to use a common magnitude and mass limit set at z = 0.45. 
Consequently, data of these four clusters have been analyzed, 
and field and group data come from these four fields in the re- 
stricted redshift range. In Table [T]useful values of the four clus- 
ters are liste d. Veloci ty dispersions (cr) were calculated using 
ROSTAT (Be ers et al.| [l990). We adopted these values to deter- 
mine cluster membership (+3cr from the cluster redshift). 

Since the projected density of cluster/supercluster members 
is low at the large clustercentric distances probed by the ICBS, 
our sample necessarily includes ~1000 "field" galaxies at red- 
shift 0.2 < z < 0.8 per survey field. This gives us an opportunity 
to compare galaxy evolution in clusters with the field over this 
epoch. 

Figure [1] shows the whole redshift distribution of the four 
fields analyzed RCS 1102, RCS 0221, SDSS 1500, and SDSS 
0845 in the redshift range considered. The cluster regions (+3cr 
from cluster redshift) are also indicated. 

Absolute magn itudes were determined using INTERREST 
dTavlor et alj|2009l) from the observed photometry. The tool, in- 
terpolating rest-frame colors from observed photometry (for de- 
tails, see Rudnick et al.,,2003.) , relates the flux in the rest-frame 
band to that through the two observed bands which bracket the 
rest-frame band. It uses a number of template spectra as guides 
and interpolates between points to determine the color relations. 
Given the apparent magnitude in the rest-frame filter, it is possi- 
ble to determine the rest-frame color. 

When photometry is available, we determine the galaxy stel- 
lar mass using the relation between M/Lg and rest-frame (B - V) 
color and the equation given in B ell & de Jong, (200 L) 

log,„(M/LB) = flfi + beiB - V). (1) 

For the B ruzual & Chario t model with an Initial Mass Function 
(IMF) of lSalpeterf ( [1951 (0.1-125 Mq) and solar metalHcity, 
Ob = -0.51 and bg - 1.45. Our broadband photometry does 
not cover the entire field of our redshift survey. If photometry 
was unavailable for a galaxy, synthetic colors were calculated 
from the flux-calibrated IMACS spectra. 

The error in the measured masses is ~0.3 dex. All our masses 
are scaled to a lKroupal (120011) IMF, by adding -0.19 dex to the 
logarithmic value of the Salpeter masses. 

The magnitude completeness limit of the ICBS is r ~ 22.5. 
At our highest ICBS redshift z ~ 0.45, we determine the value 
of the mass of a galaxy with an absolute B magnitude corre- 
sponding to r = 22.5, and a rest-frame color (B - V) ~ I, which 
is the reddest color of galaxies in ICBS clusters. In this way. 
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cluster name 


z 


cr 


^200 


^ gals 


^gals above 






(kms-') 


Mpc 




the mass lim 


SDSS0845A 


0.3308 


974.78+52.76 


2.03 


181 


100 


RCS1102B 


0.3857 


694.96+32.91 


1.40 


208 


96 


SDSS1500A 


0.4191 


527.52+37.41 


1.04 


81 


50 


RCS0221A 


0.4317 


797.63+42.50 


1.57 


201 


111 



Table 1. List of ICBS clusters analyzed in this paper, with cluster name, redshift, velocity dispersion, /?20() and number of cluster 
member galaxies (+3cr from cluster redshift). 






0.45 



Fig. 1. ICBS: redshift distribution in the four fields observed by the survey: RCS0221, RCSl 102, SDSS0845, and SDSS1500. The 
cluster regions (+3cr from cluster redshift) are also indicated (vertical dotted lines). 



the ICBS mass completeness limit at the redshifts of interest is 
M. = 1O'"-^M0. 

The color distributions of both ICBS and EDisCS samples 
(see also Figure |6] and Figure |7]i confirm that they are complete 
and unbiased down to our mass limit. 

In this paper, galaxies are given weights proportional to 
the inverse of the spectroscopic incompleteness. Since the main 
galaxy property that we wish to analyze in this work is galaxy 
stellar mass, we compute the incompleteness correction taking 
into account the number of galaxies for which there is an esti- 
mate of the mass. 

Above our mass limit, of all galaxies brighter than r = 22.5 
with spectroscopy, 91% have a mass estimates^ 

The incompleteness we have to correct for depends on both 
the magnitude and the position in the field, so it was computed 



^ 9% of galaxies have no mass estimates due to bad photometry or 
colors. 



based on the apparent magnitude and the position of each galaxy 
in the field. We subdivided each field into three different regions 
according to their distance from the centre of the main cluster of 
each field {RIRioo < 1, 1 < RIRim < 2,R/R2m > 2f| and we 
then determined the completeness weights in bins of 0.4 r mag 
around each galaxy as the ratio of the number of galaxies with 
a spectroscopic redshift and a mass estimate to the number of 
galaxies in the original photometric catalog. 



^ R200 is defined as the radius delimiting a sphere with interior mean 
density 200 times the critical density of the Universe at that redshift, 
and is commonly used as an approximation of the cluster virial radius. 
The R200 values for our structures are computed from the velocity dis- 
persions using the formula 



1.73- 



1 



(Mpc) 
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For each cluster, we excluded the Brightest Cluster Galaxy 
(BCG), identified as the most luminous galaxy, because its char- 
acteristics could alter the general conclusions. 

The final mass-limited ICRS sample of galaxies belonging 
to any environment with > lO'^^M© consists of 596 galax- 
ies. Considering also the completeness weights, the number of 
galaxies is 1295. 

3.2. EDisCS 

The multiwavelength ph otometric and spec troscopic survey of 
galaxies called EDisCS dWhite et al.l l2005h was developed to 
characterize both the clusters themselves and the galaxies within 
them. It observed 20 fields containing galaxy clusters at 0.4 < 
z< 1. 

Clusters were drawn f rom the Las Campana s Distant Cluster 
Survey (LCDCS) catalog dGonzalez et al.ll200lh . 

For all 20 fields, EDisCS consists of deep optical m ulti- 
band photometry with F0RS2/VLT jWhite et alj |2005|) and 
near-IR photometry with SOFI/NTT. ACS/HST mosaic imag- 
ing in F814W of 10 of the highest redshift clusters was 
also acquired dDesai et alj [20071) . Deep spec troscopy with 
FORS2A^LT was obtained for 18 of the fields dHallidav et alJ 
l2004tlMilvang -Jensen et al ] l2008l) . 

The FORS2 field covers the 7?2oo of a ll clusters, except fo r 
cl 7252.5-7250 where it reaches 0.57^200 (.Poggianti et al.ll2006h . 
The 7?2()o values o f our structures were co mputed from the veloc- 
ity dispersions bv'Pog gianti et al.l d2008l) . 

Photometric redshifts were computed for each object using 
both optical and infrared imaging data of the EDisCS fields 
using two independent codes: a mo dified version of the pub- 
licly available Hy perz code ( Bolzonel la. Miralles. & Pelldl2()00h 
and the code of Rudnick et al. (200lV with the modifications 
presented in Rudnick etal. (2003) and Rudnick etal. (2009). 
The accuracy of both methods is cr(&) ~ 0.05 - 0.06, where 



(km s ' ) 



RlQO 

Mpc 



6z 



zspec-zphot 



Photo-z membership (see also iDe Lucia et al.l 



I2OO4I and IDe Lucia et al. I I2OO7I for details) was established us- 
ing a modified versio n of the technique first developed in 
iBrunner & LubinI d200Q) . in which the probability of a galaxy 
to be at redshift z (P(z)) is integrated in a slice Az - +0. 1 around 
the cluster redshift to give Pdust for the two codes. A galaxy was 
rejected from the membership list if Pdust was smaller than a 
certain probability P,hresh for either code. The Pihresh value for 
each cluster was calibrated from our spectroscopic redshifts and 
was chosen to maximize the efficiency with which we can reject 
spectroscopic non-members, while retaining at least ~ 90% of 
the confirmed cluster members independent of their rest-frame 
(B-V) color or observed (V-I) color. 

For EDisCS gala xies, w e used stellar masses estimated fol- 
lowing the method of lBell & de Jong (2001) and then converted 
masses to correspond to a Kroupa (2001) IMF. Total absolute 
magnitudes were derived from photo-z fitting (Pello et al. 2009) . 
rest-frame luminosi ties were derived using the methods of 
Rudnick et al.' ('2003') and' Rudnick et all d2006h and presented in 
Rudnick et al. (2009). Stellar masses for spec troscopic members 
were also estimated using the kcorrect tool dBlanton & RoweisI 
l2007l) Fl whose masses agree with those used in this paper. For 
a detailed discussion of our mass esti mates and the consist ency 
between different methods, we refer to'Vu lcani etal] d2011'). 

For the EDisCS mass-limited sample, we used all photo- 
z members of all clusters and groups. This choice of using 
the photo-z membership instead of spectroscopically confirmed 



Clusters 
Cl 1232.5-1250 
Cl 1216.8-1201 
Cl 1138.2-1133 
Cl 1411.1-1148 
Cl 1301.7-1139 
Cl 1353.0-1137 
Cl 1354.2-1230 
Cl 1054.4-1146 
Cl 1227.9-1138 
Cl 1202.7-1224 
Cl 1059.2-1253 
Cl 1054.7-1245 
Cl 1018.8-1211 
Cl 1040.7-1155 



0.5414 
0.7943 
0.4796 
0.5195 
0.4828 
0.5882 
0.7620 
0.6972 
0.6357 
0.4240 
0.4564 
0.7498 
0.4734 
0.7043 



1080+^'" 
1018™ 

732 

' -133 

687+^^ 
666^]^^ 
648-1 

574+^2 

510!^^ 
504!"' 
486«' 
418" 



1.99 
1.61 
1.41 
1.32 
1.30 
1.19 
1.08 
0.99 
1.00 
1.07 
1.00 
0.82 
0.91 
0.70 



Groups 

Cl 1037.9-1243 
Cl 1103.7-1245b 
Cl 1103.7-1245a 
Cl 1420.3-1236 
Cl 1119.3-1129 



0.5783 
0.7031 
0.6261 
0.4962 
0.5500 



319!| 
252+^^ 
336+^'' 

-'-'"-40 

218 
166 



Table 2. List of EDisCS clusters and groups analyzed in this pa- 
per, with cluster na me, redshift, velocity d i spersion and (only for 
clusters) 7?2oo (from Hallidav et ani2004t iMilvang- Jensen et al.l 
1008; Poggi anti et alj|2008l) . 



http://cosmo.nyu.edu/mbl44/kcorrect/ 



members was made to prevent the number of galaxies being too 
low, and to permit a statistically meaningful analysis. Moreover, 
the spectroscopic magnitude limit was between 1=22 and 1=23 
depending on redshift, and the correspondin g spectroscopic stel - 
lar mass limit would be M, = 10'" *'/Wq (1 Vulcani et al.ll2010l) . 
The photo-z technique allows us to push the mass limit to much 
lower values than the spectroscopy. The magnitude complete- 
ness limit of the EDisCS photometry is 7 ~ 24 (though the com- 
pleteness remains very high to magnitudes significantly fainter 
than 7 = 24, being very close to 90% at 7 ~ 25, White et al. 
2005). We use the most distant cluster cl 1216.8-1201, z ~ 0.8 , 
to determine the value of the mass of a galaxy with an absolute 
B magnitude corresponding to 7 = 24, and a color {B-V) ~ 0.9, 
which is the reddest color for galaxies in this cluster. 

The EDisCS mass completeness limit based on photo-z's is 
then Mt - lO'^'^/W©. This is the mass limit we adopt for our 
analysis. 

BCGs were excluded from our analysis. Table|2]presents the 
list of clusters used and some basic parameters. 

The final mass-limited EDisCS sample of galaxies with 
> 10"'-2/Wo consists of 2962 objects. 



4. Definition of the environments 

We identify and separately study the different environments, 
considering the mass-limited samples of galaxies defined above. 
In the ICBS, first of all, we characterize the clusters, selecting 
only those galaxies that can be considered members, lying within 
3cr of the velocity dispersion of the cluster We then subdivide 
cluster member galaxies according to their clustercentric dis- 
tance, identifying the "cluster virialized regions" (7?/7?2(K) ^ 1) 
and the ''cluster outskirts" {R/R200 > 1). In the following, we 
sometimes subdivide the core further into three zones: inner 
parts (R/R200 < 0.2), intermediate parts (0.2< 7?/7;2oo < 0.6), 
and outer parts (0.6 < R/Rim < !)• 
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A group catalog was constructed using the standard method 
of Huchra & Geller (1982). We identify groups by a friends- 
of-friends technique, where the linking velocity distance used 
to connect friends is constant at 350 kms^\ and the projected 
linking length, Dl, scales with the incompleteness of the data as 



D, = Dn 



XM,„,„ nMu„, 
<^{M)dMI (b{M)dM 
CO %J —oo 



-1/2 



where (1) Do = OAQMpc is the linking length at a fiducial red- 
shift; (2) Zfid - 0.30, 1(r,a,6) is the incompleteness of the data 
set at a given r magnitude and position in the field, as described 
in Oemler et al. (2012a, submitted); (3) the numerator is the in- 
tegral of the galaxy luminosity function to the limiting absolute 
magnitude at the distance of the galaxy pair, corrected for galaxy 
evolution as described in Oemler et al. (2012b, submitted); and 
(4) the denominator is the integral of the galaxy luminosity func- 
tion to the absolute magnitude limit at the fiducial redshift. 

Although this method makes efficient use of all the data, it 
produces groups whose properties vary systematically with red- 
shift, because of the definition of Dl. However, we only use a 
fairly narrow redshift slice hence this drawback does not affect 
the analysis. 

We call "field galaxies" all those galaxies that are no cluster 
members. 

For EDisCS, we can identify only two main environments: 
the clusters (virialized regions and outskirts and the inner, in- 
termediate and outer regions defined as before) and the groups. 
EDisCS clusters are defined as systems with velocity dispersion 
cr > 400km i ' and their members are defined with the photo- 
z technique. Groups are defined as systems with at least eight 
spectroscopic members and velocity dispersions of I50km < 
(T < 4-OOkms^^. As for clusters, their members are defined 
with the photo-z technique. Table [3] summarizes the definitions 
adopted to characterize the different environments, separately 
for ICBS (upper panel) and EDisCS (bottom panel). 

We recall that the definitions of groups are different in the 
two samples, although, since we never directly compare them, 
this does not affect our findings. 

As akeady said, when we consider the cluster regions, in 
both samples we exclude the BCGs because they could alter the 
general trends. The properties of BCGs are in many aspects very 
different from those of other galaxies, and they are the subject 
of many studies de dicated only to this class of objects (see e.g. 
lFasanoetalJl2010h . 

In the ICBS, above its mass Hmit M, > 10'"^, on the whole, 
there are 178 cluster galaxies that belong to the virialized region, 
177 that are cluster outskirts galaxies, and 241 field galaxies, 90 
of which belong to groups (see in Table |4] the weighted num- 
bers). In EDisCS, above its mass hmit M, > 10'"-2, there are 
1268 galaxies in the cluster virialized regions, 749 galaxies in 
the outskirts, and 620 group galaxies (see Table|4]l. 

5. RESULTS 

Above the mass completeness limit, we build histograms to char- 
acterize the mass distribution of galaxies located in different 
environments. In each mass bin, we sum all galaxies belong- 
ing to the environment under consideration to obtain the total 
number of galaxies. The width of each mass bin is 0.2 dex. In 
building histograms for the ICBS, each galaxy is weighted by 
its incompleteness correction^ Errorbars on the x - axis rep- 



ICBS 


Nobs 


Nueigh! 


cluster virialized regions 


178 


339 


cluster outskirts 


177 


374 


groups 


90 


199 


field 


241 


581 


EDisCS 




^°''"m.>10I0-5Mo 


cluster virialized regions 


1268 


842 


cluster outskirts 


749 


484 


groups 


620 


421 



Table 4. Number of galaxies in each environment, above the 
completeness limit. Nobs refers to the observed number of galax- 
ies, while Nweight refers to the incompleteness-corrected num- 
bers, obtained taking into account the completeness weights. 
Upper panel: ICBS data, lower panel: EDisCS data. 



resent the width of the bin, errorbars on the y - axis are com- 
puted as a combination of the uncertainties given by the poisso- 
nian errors (Gehrels 1986) and of the uncertainties given by the 
cosmic variance^ To normalize the mass functions, we adopt a 
unique method throughout the paper and we use the total inte- 
grated stellar mass ab ove the mass completeness limit, following 
Rud nick et al ] (l20()9h . In this way, the total galaxy stellar mass 
in each histogram is equal to 1 . We adopt a normalization of this 
kind because we focus our analysis on the shape of the mass 
functions and not on the number density, which is very different 
across the different environments. 

To quantify the differences between different mass functions, 
we perform a Kolmogorov-Smirnov (K-S) test, which tells us 
whether we can disprove the null hypothesis that two data sets 
are drawn from the same parent distribution. Since the standard 
K-S test does not consider completeness when compiling the cu- 
mulative distribution (since it assigns to each object a weight 
equal to 1), we modified it, such that the relative importance of 
each galaxy in the cumulative distribution depends on its weight. 
When all galaxies have weight=l (as in the case of photo-zs) us- 
ing the modified test is equivalent to using the standard one. In 
the following, when we consider ICBS data, we always use this 
modified K-S test. 

We recall that a "positive" (statistically significant) K-S re- 
sult provides robust proof that the two distributions are different, 
but a negative K-S result does not mean that the distributions are 
identical. It is therefore useful to inspect the mass distributions 
and their upper mass, beyond the K-S test. 

In the analysis that follows, for the ICBS, we always use 
the mass limit of M. > IO'^-^Mq, while for EDisCS we use its 
proper mass limit that is M» > lO'^'^M© and sometimes also do 
some analyses using the ICBS mass cut, to qualitatively compare 
the results. 



5.1. The mass function in different environments is very 
similar 

First of all, we wish to characterize the galaxy stellar mass dis- 
tribution of galaxies located in different global environments, to 
see whether it depends on the region in which they reside. 

To begin, we compare only galaxies in the cluster virialized 



regions {R/R 



200 



1) to field and group galaxies. Using only 



^ We do not have to correct for the incompleteness of the EDisCS 
data, since we are using photo-z. 



* For field samples, we computed the errors given by cosmic vari- 
ance by estimating the differences between the field mass function and 
the mass functions obtained using separately each field. Similarly, for 
clusters we computed the errors due to cluster to cluster variance by 
comparing the mass functions of each cluster separately. 
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ICBS 



definition 



cluster virialized regions 

cluster inner part 

cluster intermediate part 

cluster outer part 
cluster outskirts 
groups 
field 



within 3o-rf„.„„, i?/i?200 < 1 
within 3o-rf„„„, RIR2Q0 < 0.2 
within 3o-rf„„„, 0.2 < R/R200 < 0.6 
within 3(rci„s,er, 0.6 < R/R200 < 1 
within 30-,,,, „f,., R/R200 > 1 
group finding Geller-Huchra method 
no cluster member galaxies 



EDisCS 



definition 



cluster virialized regions 

cluster inner part 

cluster intermediate part 

cluster outer part 
cluster outskirts 
groups 



(^struct > 400km s , R/R200 ^ 1, photo-z membership 

c.stmct > 400^m i'"', R/R200 ^ 0.2, photo-z membership 

-1 



CTstrua > 400km s 0.6 < R/R200 
o'ltmct > 400km i"' 



< 1 , photo-z membership 
R/R200 > 1 , photo-z membership 



150fcm i < (Tsiruci ^ 400km s , photo-z membership 



Table 3. Definitions adopted to characterize the several environments, for ICBS (upper panel) and EDisCS (bottom panel). 



0.3<z<0.45 



0.4<z<0.8 



,1-13 - 



'V ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I 

XICBS cluster regions (178, 339) 
~-^,*ICBS field (241, 568) 




,1-13 



' ' ' t 

X EDisCS cluster regions (1268) 
• EDisCS groups (620) 



11 11.5 
Log,„ M [M^] 




11 11.5 
Log,„ M [M^] 



Fig. 2. Observed mass function and lSchechted ( [19761) fit of galaxies in the different environments. Left panel: ICBS cluster regions 
(black crosses and solid line) and field (blue empty stars and dotted line). Right panel: EDisCS clusters (black crosses and solid 
line) and groups (red filled hexagons and dotted line). Mass functions are normalized using the total integrated stellar mass, above 
the mass completeness limit. Errorbars on the x - axis represent the width of the bin, errorbars on the y - axis are computed as 
a combination of the poissonian errors (lGehrelslll986l) and the uncertainties given by the cosmic variance and cluster to cluster 
variance. Shaded area represent Icr errors on the Schechter fits. In the labels, for ICBS both the observed and weighted numbers 
of galaxies are given, while for EDisCS only observed numbers are given. The K-S probabilities are also given. At the bottom left 
of each panel, in the M* and a space, parameters of the Schechter fit are shown. Errorbars represent the Icr errors. No significant 
differences are evident between the galaxy stellar mass distributions in different environments. For ICBS, similar results are also 
obtained comparing the cluster regions and the non-cluster regions (plot not shown). 



ICBS data, we can contrast the most widely different environ- 
ments: the cluster virialized regions and the field. In the left 
panel of Figure |2] we can see that no significant differences 
are evident between the shapes of the galaxy stellar mass dis- 
tributions for the two environments. The K-S test is unable to 
reject the null hypothesis that the two distributions are drawn 
from the same sample (Pk-s ~ 46%). The mass function agrees 
well with the K-S result: its shape is similar, within the errors, 
in the two environments. To increase the quality of the data 
statistics, we then compare the cluster virialized region with the 
field+outskirts galaxies. Again, we do not find detectable differ- 
ences (Pk-s ~ 40%, plot not shown). 



In the right panel of Figure |2] using EDisCS data we can 
characterize clusters and groups, reaching even lower masses 
(M, ~ lO'^-^Mo). Also in this case, there is no evidence of a 
dependence of the shape of the mass function on the environ- 
ment. The Pk-s is inconclusive (Pk-s ~ 21%) with a high sta- 
tistical certainty because of the large number of galaxies. From a 
visual analysis of the plot, it is clear that the shapes of the mass 
functions are very similar. 

In principle, photometric red shifts might be very u ncertain, 
especially for blue galaxies ('see iRudnick et all 1200 91 and the 
errors could bias the derived stellar mass functio ns for EDisCS. 
Howev er, we tested that the magnitudes at which iRudnick et al] 
(l2009h found the discrepancy correspond to lower stellar masses 
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than those considered in the paper: our Mg is always < -20.2. In 
addition, the photo-z counts are fully consistent with the statisti- 
cally background subtracted counts for both red and blue galax- 
ies and, in the mass range in common, the mass function de- 
termined from photo-z's is in agreement within the errors with 
the mass function determined using only spectroscopic mem- 
bers and taking into account spectroscopic completeness weights 
dVulc ani et al. 201 1). The results are compatible and robust and 
this gives us confidence in the photo-z technique. 

For both ICBS and EDisCS, results in the low mass end are 
very trusty and reliable, while in the high mass end the uncer- 
tainties do not allow us to definitely characterize the trend of the 
mass functions, hence we can not state whether the similarities 
are real or simply due to the size of the errors. 

To give strength to our results, we also perform a fit of the 
mass functions using the least square fitting method. Assuming 
that the number density <i>(M) of galaxies is described by a 
ISchechte^ (119761) function, the galaxy stellar mass function then 
is 

a)(M) = (In 10) X <D* X [io(M-M')(i+a)j ^ exp[-10^'^-'^'>] (2) 

where M = log(M,/M0), a is the low-mass-end slope, M* - 
log(M* /Mq) is the characteristic stellar mass at which the mass 
function exhibits a rapid change in the slope, and O* is the nor- 
malization. Schechter functions are fit to galaxies only above our 

conservative completeness l imit. 

Table|5]gives the best-fit lSchechteii d 1 9761) parameters for the 
mass functions of galaxies in different environments (see also 
£21and jO . 

Leaving free all the parameters of the fit, we find that mass 
functions in different environments show comparable parame- 
ters, within the Icr errorQ This supports the finding that the 
shapes of the mass functions seem not to depend on the global 
environment. This is reliable above all for EDisCS, for which 
the statistics is high and hence the parameters are rather well 
constrained^ 

We stress that our comparison is really about the shape (here 
defined by logM* and a) and not about the normalization O*, 
which is arbitrary and almost never comes into play. 

We have to note that, as seen in other works, e.g., iBeif et alJ 
l2003HBaldrv et alJl2008 ) the Schechter fit is not able to properly 
describe the very high mass end of cluster and group mass func- 
tions, which both show a sort of bump at M, ~ lO'^'^M©, how- 
ever it well fits the distributions for M» ~ IO'^Mq. In EDisCS, 
since we are using photo-z memberships, the bump might be due 
to a contamination by interlopers, however, we have checked out 
that most of the photo-z members are also spectroscopic mem- 
bers and that galaxies in the bump do not show bigger photo-z 
errors than lower mass galaxies. Finally, we note that a similar 
bump has also been detected in the spectroscopic low-z clus- 
ter sample WINGS (Wide-field Nearby Galaxy-cluster Survey - 
[Fasano et al. 2006) (see S2|. 

For ICBS, we can also separately consider narrowly defined 
environments, to see whether there is any difference among clus- 

' Taking into account the fact that the two parameters are correlated, 
we explored a grid of a and M' parameters, finding the corresponding 
X~ values and the likelihood of having the same couple of values and 
found that mass functions are in agreement within Icr. 

For EDisCS, we also comput ed the best-fi t ISchechtej ( Il976h 
parameters using the STY ( Sanda ge"et all Il979h method (see, e.g, 
[Marchesini et al. 2009) which is an unbinned maximum likelihood 
method and the parameters are compatible within Icr error We do not 
adopt this method throughout the entire paper because it is not trivial to 
take into account ICBS' weights. 



ter, group, and field galaxies without groups (plot not shown). 
Even though there is moderate statistical uncertainty, and the 
mass functions are hence noisier, their shapes are not obviously 
different (the Pk-s is always inconclusive). 

The robustness of our results is demonstrated by the similar 
lack of environmental dependence of the mass functions for both 
EDisCS and ICBS0 

These results s eem not be in line with the findings of 
iKovac et"an (l2010h . who studied ~ 8500 galaxies from the 
zCOSMOS-bright redshift survey in the COSMOS field. They 
found that the stellar mass function exhibits a different shape for 
samples of galaxies in different environments (groups, field, and 
isolated galaxies) at least up to z ~ 0.7. Their stellar mass func- 
tion shows an upturn at low masses in the group environment and 
they found that more massive galaxies preferentially reside in the 
groups. However, their samples are much deeper (they go down 
to Mi, ~ IO^'^Mq) hence they can inspect lower mass galaxies 
and above all the definitions adopted to select group and field 
galaxies are not comparable to ours. 

We note that, even if the mass functions in different envi- 
ronments appear to have similar shapes, there is a hint that they 
extend up to different maximum masses (the so called mass func- 
tion cut-off): for the ICBS, the most massive galaxies in clus- 
ters have M, ~ lO^ '^Mo, in the field M. ~ 10"-^Mo, and in 
groups M» ~ 10"^Mo. In EDisCS, clusters virialized regions 
can contain galaxies as massive as M, ~ IO'^'^'Mq, cluster out- 
skirts M. ~ 10'2 'Mo and groups M. ~ IO'^-^Mq. 

The median masses are not significantly different between 
the diff'erent environments: in the ICBS, clusters, groups, and 
field all have a value close to M, ~ 10'"^Mo above the ICBS 
limit, while EDisCS clusters and groups have median masses 
very close to M, ~ 10'° ''^ Mq above the EDisCS limit. 

5.2. The evolution of the mass functions is very similar in 
different environments 

As we have seen in the previous section, at least for galaxies 
with Mt > IO'^'^^'^'^Mq the shape of the galaxy stellar mass 
function does not seem to depend on the global environment in 
which galaxies reside. 

In Calvi et al. (in preparation), we have carried out a similar 
analysis of a mass-limited sample (M» > 10'"'^^Mq) of galaxies 
in the local Universe, using the Padova Millennium Galaxy and 
Group Catalog (PM2GC) (ICalvi. Poggianti. & Vulca"nil 1201 ih 
and WINGS. As in this work, in the local Universe we have 
found that, excluding the BCGs, clusters, groups, and the field at 
low-z follow comparable mass functionsl3 

Having found similar results at both redshifts, we wish to in- 
vestigate whether the evolution of the galaxy stellar mass func- 
tion c hanges with e nvironment. 

In I Vulcani et al.l (1201 II) . we compared cluster galaxy stellar 
mass functions at low and high-z using WINGS and EDisCS 
data and found a strong evolution from z ~ 0.8 to z ~ 0, which 
we attributed primarily to mass growth due to star formation in 
late-type cluster and infalling galaxies. We found that the shape 
at Mt > 10"Mq does not evolve but that the mass function at 
high redshift is flat below M» ~ IO^'^Mq, while in the Local 



' We recall that ICBS is a spectroscopic sample, hence the separation 
into the several environments is very reliable, EDisCS provide a large 
sample and also it extends to lower masses. 

However, at low-z the mass functions in the different environments 
span different ranges of masses, as also happens to a certain extent at 
higher z (see Calvi et al. in preparation). 
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Table 5. Best-fit lSchechteil (1 19761) parameters (a, M*, (p* ) for the mass functions of galaxies in different environments and at different 
redshift (see ^5^21 and 1531). 







logM*(logM,/Mo) 


a 


(I)*(log</>) 


ICBS 


cluster regions 


11.37±0.28 


-1.29+0.41 


30.54+27.72 




cluster outskirts 


11.40±0.11 


-1.29+0.00 


33.07+6.70 




field 


11.00±0.12 


-0.89+0.29 


134.28+40.23 


EDisCS 


cluster regions 


11.15±0.07 


-1.03+0.08 


142.53+23.70 




cluster outskirts 


11.25+0.06 


-1.03+0.00 


73.71+5.10 




groups 


11.18+0.06 


-0.97+0.08 


70.99+9.99 


WINGS 


cluster regions 


10.82+0.13 


-0.88+0.31 


219.48+74.79 


PM2GC 


general field 


10.96+0.06 


-1.12+0.12 


173.65+32.11 




11 11.5 

Log,„ M [Mj 

Fig. 3. Comparison between the mass function in clusters and in 
the field at different redshifts. Low-z: red filled squares: PM2GC 
(general field), cyan empty squares: Baldry et al. (2008) (general 
field), blue empty triangles: WINGS (clusters). Intermediate-z: 
black crosses: ICBS (field), black skeleton tria ngles: ICB S (clus- 
ters). Higher-z: magenta skeleton stars : Poz zetti et al.l (1201 Ol) 
(field), green filled triangles EDisCS (clusters). Mass functions 
are normalized using the the total integrated stellar mass above 
log M»/M0 > 11.2. Errorbars and the inset are as in Figure |2] 
The evolution of the shape of the mass function seems not to de- 
pend on environment, being similar in clusters and in the field. 
However, the amount of growth as a function of cosmic time is 
not the same at the low- and high-mass end. 



Universe it flattens out at significantly lower masses. The pop- 
ulation of M* - lO'"'^ - 10'" ** galaxies must have grown 
sig nificantly betw e en z = 0.8 and z = 0. 

iPozzetti et al] (1201 Oh . using data from zCOSMOS-bright 
10k spectroscopic sample, quantified the evolution of the mass 
function in the field, from z ~ 1. They used data from 
[Baldrv et al. (2008), who selected galaxies from the NewYork 
University Value-Added Galaxy Catalog sample, as reference at 
z = 0. They found a continuous increase with time in the mass 
function for log M/M, < 11, while a much slower increase at 
higher masses. 

We are now in the position to compare the evolution of the 
mass function in clusters with that in the field. In Figure [3] we 
examine the evolution from z ~ 0.4 - 0.8 to z ~ 0. We show the 



mass functions for WINGS clusters and lBaldrv et alj (l2008l) and 
PM2GC general field at z ~ 0.07, those of clusters and field from 
ICBS at z ~ 0.4 and those of EDisCS clusters and Pozzetti et al| 
( 1201 Oh field (private communication) at z ~ 0.6. Perhaps in con- 
trast to expectations, the evolution does not depend on global en- 
vironment, being similar in clusters and the field. The shape of 
the mass functions of the field and clusters overlap quite consid- 
erably, within the eiTors. As cosmic time goes by, the number of 
galaxies at low-to-intermediate mass grows proportionally with 
respect to the number of massive galaxies, in the same way in 
clusters and in the field. 

However, it is hard to quantify the evolution from z ~ 0.4 
and z ~ 0.6 to z - 0, du e to th e size of the errorbars. 
In Table |5] best-fit ISchechteii d 19761) are given for WINGS, 
PM2GC, ICBS and EDisCS samples, for which stellar mass es- 
timates and mass functions are derived in the same ways and 
hence results are comparable. Instead, we decided not to quan- 
titat ively compar e our findings to those of iPozzetti et alJ (1201 Oh 
and lBaldrvetal] {2008) since results obtained using heteroge- 
nous data and slightly different redshift ranges are hardly com- 
parable to draw definite conclusions. Indeed, we stress that us- 
ing stellar mass functions derived from independent works can 
be dangerous, as the different methods to derive stellar masses, 
with different assumptions, can results in differences in the mass 
estimates by a factor of ~2-3. The effect of these systematic dif- 
ferences can be significant at the high-mass end, biasing the de- 
rived evolution. 

The best-fit parameters for WINGS and PM2GC well agree 
within the errors, as EDisCS and ICBS do, even though the large 
uncertainties of the ICBS sample do not allow to well constrain 
the mass function. Going from lower to higher z, the character- 
istic mass increases, while a clear trend can not be defined for 
the low mass end slope. 



5.3. The mass function in different cluster regions 

As we considered clusters, groups, and the field separately, we 
now shift our attention to clusters alone, to analyze more care- 
fully the different regions of clusters and compare the central 
parts with the outer regions. 

As described in ^ we can subdivide clusters into several 
regions at different galaxy clustercentric distances. 

We first compare the cluster regions {R/R200 < 1) with the 
outskirts {R/R20Q > 1) (Figure |4]i. Using both the ICBS and the 
EDisCS data sets and applying the K-S test, we are unable to de- 
tect any difference (Pk-s ^ 10% in both cases). This is also sup- 
ported by the analysis of the Schechter fits (see Table|5]l. In this 
case, we choose to fix the low-mass end slope of the mass func- 
tion of galaxies in the cluster outskirts and we adopt the a value 
we found for the virialized regions. Again, the fit is not able to 
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0.3<z<0.45 0.4<z<0.8 




Log,„ M [M^] Log,„ M [M^] 

Fig. 4. Observed mass function and lSchechteii ( Il976h fit of galaxies at different clustercentric distances (R/R200 < 1 and R/R200 > 1 ) 
for ICBS (left panel) and EDisCS (right panel). Black crosses and solid lines represent cluster regions, and red empty stars and dotted 
lines the cluster outskirts. Mass function normalizations, errorbars, labels and insets are as in Figure|2] In both panels, no statistically 
meaningful differences are detected between the mass functions of galaxies located at different clustercentric distances. 



describe the bump observed in EDisCS at at M» ~ IO'^ -'Mq, 
but it works well at lower masses. The parameters of the fit sup- 
port the finding of similarity of the shape of the mass functions, 
seen in both panels in Figure 31 where despite the poor statistics 
of the ICBS sample, it is clear that the mass functions are very 
similar 

We next compare the three different zones of the virialized 
regions, i.e. the inner (0 < R/R2Q0 < 0.2), the intermediate 
(0.2 < R/R200 < 0.6), and the outer (0.6 < R/R200 < 1) parts. 
Both for ICBS data and EDisCS, the Schechter parameters are 
compatible within the errors and the K-S test is always incon- 
clusive and it is unable to detect any strong variation with clus- 
tercentric distance. Figure [5] shows the mass functions for the 
inner and outer regions for EDisCS dataH 

To summarize, no overall differences are detected between 
the cluster virialized regions and outskirts and this agrees with 
the previous result that global environment does not alter the 
shape of the mass distribution. We can consider the outskirts as 
a transition region between the cluster virialized regions and the 
field and no differen ces are detected between the m. Our results 
agree with those of Ivon der Linden et a derived using 

SDSS data: excluding the BCGs, they showed that there is no 
evidence for mass segregation in clusters, the median mass of 
cluster galaxies being invariant with cluster radius. 

5.4. The mass function of red and blue galaxies does not 
depend on environment 

In the previous subsections, we have shown that there appears 
to be no dependence of the shape of the mass function on the 
global environment. Our finding is quite surprising, because it 
is known that galaxies located in different environments have 
different distributions of morphological and star formation prop- 



" We do not show the plot for ICBS data since the quality of the 
statistics is quite low. 



0.4<z<0.8 




Log,„ M [Mj 

Fig. 5. EDisCS: Observed mass functions and' Schechted (Il976h 
fit of galaxies in dfferent cluster regions within R/R200 - 1- 
For sake of clarity, only inner (black crosses and solid line) and 
outer parts (red filled hexagons and dotted line) are plotted. Mass 
function normalizations, errorbars, labels and the inset are as in 
Figure 12] The shape of the mass function is similar in the differ- 
ent regions. 



erties, hence a natural question is then whether different galaxy 
types also follow the same mass distribution. 

We decided to subdivide galaxies by color, into red and 
blue, to separate galaxies with diff erent star forma tion proper- 
ties. Adopting the cut proposed bv iPeng et alj ( 12010 ) and con- 
verting to our adopted IMF and to the Vega system, galaxies are 
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Fig. 6. Stellar mass versus rest-frame U-B color for our whole 
mass complete sample in ICBS. The line separating red and blue 
galaxies is also plotted. 
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Fig. 7. Stellar mass versus rest-frame U-B color for our whole 
mass complete sample in EDisCS split into four redshift bins. 
The line separating red and blue galaxies is also plotted. 

assigned to the red sequence using the cut 

{U - B)vega > 1.10 + 0.075 X log( ^(p^^ ) - 0.18 x z - 0.88, 

while the other galaxies are assigned to the blue cloud. 

Figure |6]and Figure |7]show the the rest-frame (U-B) color 
as a function of stellar mass estimates and the cut adopted to sep- 
arate the red and blue populations for the samples of ICBS and 
EDisCS, respectively. Since EDisCS covers a quite wide range 
of z, we divided the whole sample into four redshift slices, to en- 
able us to more clearly visualize the red sequence and the blue 
cloud. 

At 0.4 < z < 0.8 contamination from dusty, star-forming 
galaxies can be quite important when selected red galaxies using 
only one color As a consequence, we stress that the stellar mass 
function of red galaxies is not the stellar mass function of quies- 
cent galaxies, but it is partly contaminated by dusty star-forming 



Table 8. Fractions red galaxies as a function of mass both in 
EDisCS and ICBS samples. Errors are computed as binomial 
errors. For the ICBS, only weighted fractions are given. 
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galaxies. In EDisCS, the level of contamination of the adopted 
red galaxy selection fr om dusty star-forming objects found using 
Bramm er et al.l (l20Tll) method is ~10%, while in ICBS we can 
not compute it due to the limited wavelength coverage. 

Tables |6] and [T] show the fraction of red and blue galax- 
ies in the two samples. As expected, these fractions strongly 
depend on environment. In the cluster regions, red galaxies 
dominate the whole population, above all at higher masses: for 
M, > 10"'^Mo, about 90% of all galaxies in the ICBS are red, 
and about 60% in EDisCS. EDisCS and ICBS fractions differ but 
they may not be directly comparable, in particular because they 
cover different redshift ranges. For the EDisCS mass cut, which 
is lower, in EDisCS the red fraction is slightly lower (~ 54%), in- 
dicating that low mass galaxies are preferentially blue. In ICBS 
the red fraction reaches a minimum of 41.7% in the cluster out- 
skirts. 

For both samples, we also computed the fraction of red 
galaxies in different mass bins, finding that it strongly depends 
on mass, hence indicating that these galaxies dominate the high 
mass end, in any environment. These fractions (see Table [8]) are 
consistent within \-2cr with the fractions of quiescent galax- 
ies ob tained from the number densities given by Brammer et all 
(|201I|) for galaxies in the field at z ~ 0.613: 50 ± 12% for 
log MJMe = 10.2-10.6,75±19%forlogM./Mo = 10.6-11, 
and 92 + 33% for log M./Mq = 11- 1 1.6n3 

In Figure |8] we show the mass function of red and blue 
galaxies using EDisCS data. Above M, > lO'^ ^M©, in clusters 
(left panel), in the outskirts (central panel) and in groups (right 
panel), blue and red galaxies have different mass distributions 
{Pk-s ~ 0% in all cases). This is also immediately clear when 
looking at the plots: in all cases, blue galaxies tend to have pro- 
portionally more low mass galaxies than red galaxies, especially 
in clusters (both virialized regions and outskirts). The blue mass 
function is therefore steeper than the red one. This is the only 
case where in this paper we detect a difference between the mass 
functions compared. 

Red and blue galaxies also have different mass functions in 
the ICBS field and outskirts (plots not shown). The quality of 
the statistics of the ICBS clusters instead is too low to detect 
any difference. This is clearly illustrated when we extract from 
the EDisCS sampleP^ a subsample that has the same number of 



We note that the adopted choice to subdivide passive/star-forming 
galaxies strongly influences the fractions, so every comparison be- 
t ween d ifferent cuts has to be analyzed very carefully. In addition, 
iBrammer et al.l ( 201 1 ) give the quiescent/star-forming number densities 
only for field galaxies in a redshift range slightly higher than the ICBS' 
one. 

" In this case, we use the EDisCS subsample with M, > lO'^^M©. 
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Table 6. Fractions of blue and red galaxies in the ICBS sample (0.3 < z < 0.45). Errors are computed as binomial errors. Both 
observed and completeness-weighted numbers are listed. 
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Table 7. Fractions of blue and red galaxies in the EDisCS sample (0.4 < z < 0.8). Errors are computed as binomial errors. Both the 
EDisCS proper mass limit (M, > IO'^^Mq) and the ICBS mass limit (M, > IO'^-^Mq) are considered. 
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galaxies as in the ICBS: in this case the K-S test is not able any- 
more to detect any difference between the red and blue popula- 
tion. Performing 1000 Monte Carlo simulations, this test allows 
us to conclude that the ICBS clusters result are not reliable, but 
simply biased due to the poor statistics. 

Moreover, when we adopt the ICBS mass cut for EDisCS 
cluster regions, we still find that red and blue galaxies follow 
different mass distributions (Pa:-s=0.18%). Hence, differences 
are not limited to those found at low masses. 

Finally, we can separately compare the mass function of blue 
(Figure |9]) and red galaxies (Figure [TOt in different environ- 
ments. From the K-S and the plots there is generally a lack of an 
environmental dependence for both blue galaxies and red galax- 
ies. 

We also derive the best-fit ISchechted (Il976h parameters 
(Table|9]). In this case, the data lack the number statistics to allow 
robust estimation of the faint end slope a. A number of strate- 
gies could have been adopted; we choose to assume a fixed value 
of a. Fo r red galaxies, we assume a = -0.7 (the same value 
found bv lBorch et al.ll2006h . while for blue galax ies we assume 
a = -1.25, (the value used bv lBorch et aDdioOd) (alphas- 1.45) 
is too high to well describe our mass functions). Our choice is 
also supported by the fact that for EDisCS, when we leave free 
all the parameters, we recover very similar values for a. Due 
to both the high mass completeness limit and the low number 
statistics, we are not able to find the best-fit parameters to de- 
scribe blue galaxies in the ICBS clusters. Again, we note that 
the Schechter fit does not describe the very high mass end of 
EDisCS galaxies. 

The analysis of the parameters supports the evidence that in 
all environments, red and blue galaxies, while they have different 
slopes at low masses (a values are different), they show compa- 
rable values of M*, within l-2cr errors, indicating that actually 
the shape of the high mass end is rather similar for all galaxies. 
Our findings are in agreement with the analysis performed by 
[Borch et al. (2006) and Ilbert et al. (2010), who also found that 
at intermediate redshifts galaxies of different color show rather 
similar M*[3 



At z ~ 0.5, iBorch et al.l ( 120061) found that red g alaxies have M' = 
jQio.95±o.io gjjjj galaxies have M* = 10'°''^^'"^. lirbert et alJ (l2010l) 



To conclude, the analysis of the best-fit parameters confirms 
that galaxies of the same color are described by mass functions 
with similar shape in all the environments: the same a can be 
chosen and consequently the determined M*s are compatible in 
all environments at l-2cr level. The results are particularly robust 
for EDisCS, given the high number statistics. 

6. Discussion 

Results coming from EDisCS and ICBS samples are generally in 
agreement, or at least compatible. Unfortunately, we only quali- 
tatively compared our findings, but we couldn't perform any di- 
rect comparison, because the two surveys have different charac- 
teristics, their data not being homogeneous and covering slightly 
different redshift ranges. 

By sampling the mass functions for such different global 
environments, we study dark matter haloes of a wide range of 
masses, and yet, find that there is no obvious dependence of 
the galaxy mass distribution. In general, the results of this work 
suggest that the galaxy stellar mass function do not depend on 
the global environment: galaxies located in clusters, groups, and 
field seem to follow mass distributions characterized by similar 
shapes. This result is surprising and, at some level, perhaps con- 
trary to most expectations. 

We will need to compare our results with theoretical ex- 
pectations, to understand whether simulations predict any mass 
segregation with environment, considering both the initial and 
evolved halo mass and how they predict the evolution with red- 
shift as a function of the environment (Vulcani et al. in prepara- 
tion). 

We need now to discuss several aspects, in order to clarify 
the emerging picture. 

6.1. The evolution of the mass function in different 
environments 

In lVulcani etll] (|20TTI) . we argued that the evolution observed 
in clusters is driven by the mass growth of galaxies caused by 



found that red sequence galaxies have M' = lo"'"*'"'^ and intermedi- 
ate activity galaxies have M* = lo'^'^ioos 



13 



Benedetta Vulcani et al.: No dependance of the mass function on global environment 



0.4<z<0.8 0.4<z<0.8 0.4<z<0.8 




10 10.5 11 11.5 12 12.5 10 10.5 11 11.5 12 12.5 10 10.5 11 11.5 12 12.5 



Fig. 8. EDisCS: Observed mass function and ISchechteii (1 1 9761) fit of blue and red galaxies in clusters (left panel), in the outskirts 
(central panel) and in groups (right panel). Blue filled squares and solid lines represent blue galaxies, and red empty triangles 
and dotted lines red galaxies. Mass function normalizations, errorbars, labels and insets are as in Figure |2l In the insets, errorbars 
represent the Icr errors (solid line) and 2cr errors (dotted line). In clusters, in the outskirts, and in groups, blue and red galaxies have 
different mass distributions. 



0.3<z<0.45 0.4<z<0.8 




Log,„ M [M^] Log,„ M [M^] 

Fig. 9. Observed mass function and ISchechterl d 19761) fit of blue galaxies in the different environments in ICBS (left panel) and 
EDisCS (right panel) samples. Black crosses and solid lines represent the cluster regions, blue filled squares and dotted lines the 
cluster outskirts, and green empty triangles and dashed lines the groups. For ICBS cluster blue galaxies not fit has been found due 
to the small number of the bins. Mass function normalizations, errorbars, labels and insets are as in Figure|2l In the insets, errorbars 
represent the Icr errors (solid line) and 2cr errors (dotted line). For blue galaxies, no differences can be detected between the mass 
functions of galaxies located in clusters, groups, and in the field. 



star formation in both cluster galaxies and, most of all, in galax- 
ies infalling from the cluster surrounding areas. In that prelim- 
inary analysis performed using inhomogeneous data, we also 
hypothesized that infalling galaxies could in principle follow a 
different (steeper) mass distribution (environmental mass segre- 
gation) than cluster galaxies, and hence give a major contribu- 
tion to the intermediate-to-low population, although we found 
no evidence of any difference between the cluster mass func- 
tion and field mass functions taken from the literature. In this 
work, we have been able to analyze the mass function of galax- 
ies in the field and, most importantly, in the cluster surrounding 
areas. As a consequence, we can characterize the mass distribu- 
tion of galaxies that are supposed to fall into clusters between 



high and low redshift and compare it to the cluster mass func- 
tion at the same redshift. In summary, we have shown that, at 
least for the mass ranges considered (log M^/M© > 10.5), the 
mass function is invariant with the environment and that galaxies 
located in different environments follow very similar mass dis- 
tributions. Hence, the observed evolution of the mass function 
in clusters (I Vulcani et al.ll20lTI) probably can not be explained 
by galaxies of different masses residing in independent/different 
environments, at least when we employ a quite high-mass cut. 
Star formation is the remaining major process left to explain the 
observed mass growth, both in clusters and in the field. 

Moreover, by analyzing also the field mass function of galax- 
ies in the local Universe (Calvi et al. in preparation), we have 
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Fig, 10. Observed mass function and ISchechteii (fT976l) fit of red galaxies in the different environments in ICBS (left panel) and 
EDisCS (right panel) samples. Black crosses and solid lines represent the cluster regions, red filled squares and dotted lines the 
cluster outskirts, and green empty triangles and dashed lines the groups. Mass function normalizations, errorbars, labels and insets 
are as in Figure |2l In the insets, errorbars represent the Icr errors (solid line) and 2cr errors (dotted line). For red galaxies, no 
differences can be detected between the mass functions of galaxies located in clusters, groups, and in the field. 

Table 9. Best-fit [Schechted (1 19761) parameters (a, M*, <p*) for the mass functions of galaxies in different environments and of 
different colours. 
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investigated the evolution of the mass function in the field from 
redshift z ~ 0.4 to z ~ and com pared this to the evolution found 
in clusters ( Vulcani et al.l l201 1). Our results show that, at least 
for log M,/M0 > 10.2, the galaxy stellar mass function evolves 
in the same way in all environments, but not of a fixed growth 
at both ends: as time goes by, it becomes steeper, indicating that 
the number of intermediate-mass galaxies grows proportionally, 
in the same way in the field and in clusters. 

The evolution of the mass function of cluster galaxies might 
also be driven by the change in the mass function of the infalling 
galaxies that make the clusters (hence by the field mass func- 
tion). Nevertheless, our finding is quite surprising, because it 
is well known that galaxies in different environments and with 
different stellar masses have different star formation properties 
and are subject to different mechanisms. In clusters and in the 
field we expected that the processes that suppress or halt star 
formation were different, and hence that the mass growth was 
different in different environments and had different time scales. 
At these redshifts, most of the galaxy mass instead appears to 



have already been assembled and that environment-dependent 
processes have had no significant influence on galaxy mass. The 
star formation-mass relation has a similar trend in clusters and 
in the field, although a lower median SFR (by a factor of ~ 1 .5) 
is detected for clust er star-forming galaxies than for the field 
(I Vulcani et al.ll2010 l). In any event, the slightly different mass 
growth in the different environments appears to be insufficient to 
considerably alter the mass distributions. 

6.2. The blue and red mass functions 

In ^5.41 we focused our attention on red and blue galaxies, to 
see whether the shapes of the mass function of galaxies with 
different star formation properties varies with the environment. 
In summary, we have generally found that in each environment 
red and blue galaxies are regulated by different mass functions 
(where that of blue galaxies is always steeper than that of red 
galaxies - see Figure [8]l, while blue and red galaxies separately 
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follow almost the same mass function in all environments (see 
Figure |9]and Figure [TOl l. 

However, we have found that the fraction of red and blue 
galaxies strongly depends on environment: in the cluster regions, 
red galaxies dominate the galaxy population, while blue galaxies 
are found mostly in the field. Therefore, it is quite surprising that 
the total mass function is almost always the same in all environ- 
ments. 

Figure [TT] shows the fit of the mass functions (see Tables |5] 
and|9l) for all, red and blue galaxies in the different environments, 
extrapolated toward lower masses, without applying any normal- 
ization. The Figure shows well how blue galaxies dominate in 
number the mass functions at low masses, in every environment. 
Thus, the shape of total mass function at low masses is regulated 
by the shape of the blue mass functions, in all environments. 
From the analysis of the fits, we have found that blue galaxies 
can be described by adopting the same value of a, and this ex- 
plains the similarity of the total mass functions in the different 
environments, for masses < M*. In addition, the similarity of M* 
for red and blue galaxies in clusters, groups and field, makes the 
high mass end of the mass functions similar in all environments. 

To conclude, the lack of dependence of the total mass func- 
tion on environments can be described with a careful analysis 
of the Schechter parameters used to describe the red and blue 
galaxies in the different environments. 

6.3. Global and local environment 



In IVulcani et alj (1201 2l) . we have specifically analyzed the de- 
pendence of the mass function on local density, using collec- 
tively the WINGS, PM2GC, ICBS, and EDisCS samples. We 
have found that in both the local and distant Universe, in both 
clusters, and the general field, local density plays an important 
role in driving the mass distribution: in general, lower density 
regions host proportionally a larger number of low-mass galax- 
ies than higher density regions. In particular, local density in the 
general field regulates the shape of the mass function at both 
low and high masses. In contrast, the situation in clusters is 
slightly different: local density is important only when we can 
reach a quite low-mass cut (logM^/Mg < 10.1 in WINGS and 
log M./Mq < 10.4 in EDisCS), while for high-mass thresh- 
olds we detect no influence. We also found that not only the 
shape of the mass function is different, but also the highest mass 
reached: very massive galaxies are located only in very high den- 
sity regions, while they are absent in the the lowest density re- 
gion studied (the so-called mass segregation). Hence, above the 
same mass, we detected differences among mass distributions 
of galaxies located at different local densities but not of those 
of galaxies in different global environments. To summarize, if 
we put together our results, they suggest that global and local 
environment seem to have a different influence in shaping the 
mass functions. While the global environment seems to be less 
influential, local density is important in determining the most 
fundamental of all galaxy properties, the galaxy mass. Global 
and local environment are clearly related to different physical 
processes, and their different behavior in influencing the mass 
function is important to understanding the drivers of galaxy for- 
mation and evolution. 

6.4. Some caveats 

We reemphasize that our results are valid only for the mass 
ranges covered by our samples: we do not know anything beyond 



log M,/Mo < 10.2 in clusters and groups and log M./Mq < 10.5 
for the field. In principle it could be possible that at lower masses 
the situation is very different. Surveys with lower mass com- 
pleteness limits will be needed to establish the role of the envi- 
ronment for low-mass galaxies. In all cases, we have to keep in 
mind that the mass limit that we have adopted is not so low: the 
mass limit of EDisCS, of M, ~ 1.6 x 10"^Mo, is beyond the mass 
limit proposed by Kauffmann et al. (2003) (3 x IO'^Mq) to sub- 
divide low-z galaxies characterized by different properties, such 
as age of the stellar populations, surface mass density, concen- 
tration, and star formation rates. Moreover, above the same mass 
limit we observe variations in the mass function both depend- 
ing on galaxy colours and in different local environments (see 
IVulcani et al. l l20Tl . indicating again that the mass Umit is not a 
limiting factor. 

Moreover, the lack of evidence for a dependence of the mass 
function on global environment could also be partly due to small 
number statistics, especially in some cases. Most of our results 
are reliable, above all those that are confirmed by both our sam- 
ples (the statistical certainty is quite high for all plots of the 
EDisCS sample), but it could be very important to have, for 
example, a larger spectroscopic sample of field galaxies, to as- 
sess whether isolated galaxies still follow the same mass dis- 
tribution. In addition, we have to remember that our EDisCS 
sample has been selected using photo-z data, hence the level of 
contamination, even if very low (see e.g. Hallidav et al. 200^ 
iMilvang- Jensen et al.ll2008l) . may not be totally negligible. As a 
consequence, studies exploiting the capabilities of larger spec- 
troscopic surveys are needed to confirm our results. 

Our results also depend on our adopted IMF: we implic- 
itly assume that it is universal, regardless of time, environment, 
and galaxy morphological types. Of cou rse, this may not al- 
ways be the case; for example, recently, iGunawardhana et all 
(1201 ih argued that there could be a dependence of the IMF on 
the star formation rate: galaxies with a high absolute value of 
their SFR may have an IMFs with flatter power-law slopes than 
galaxies with low star formation rates. Recent results have also 
found that the IMF seems to depend on the velocity dispersion 
of the galaxy, with a bottom-heavy IMF for high-velocity dis 
persion objects, and a Kroupa-l ike IMF for low-v elocity dis 
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persio n objects (e.g . . iTreu et al.ll2009: .van Dok kum & Conrov 



Dutton et al.ll201 it ICappell ari et al 



201 1 : iThomas et al 

i2012[lSpiniello et al.ll2012h . If this is the case and the IMF is not 
universal, the results could be quite different. We could expect 
that, if the IMF varies as proposed, the mass functions would be 
"stretched": the low mass part would largely remain unchanged, 
while very massive galaxies would be even more massive. 

However, it is difficult to envisage how our results could be 
due to a conspiracy of the IMF, which should be different for 
galaxies in different environments and of different morphologi- 
cal types in a way that the total mass function in clusters, groups, 
and field is the same. 



7. Conclusions 

In this paper we have analyzed the shape of the stellar galaxy 
mass functions in different environments (mainly clusters, 
groups, and field). We have studied two different ranges of in- 
termediate redshifts: using ICBS we considered 0.3 < z < 0.45, 
and using EDisCS we considered 0.4 < z < 0.8. We have stud- 
ied mass-limited samples with M, > IO'^-^Mq for ICBS and 
M, > 1O"*^M0 for EDisCS, hence we could characterize only 
the relatively massive end of the mass functions. 
Our main results can be summarized as follows: 
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Fig. 11. ISchechteJ (Il976h fit of blue (blue dotted lines) and red (red dashed lines) galaxies for ICBS (upper panels) and EDisCS 
(lower panels) in clusters (left panels), in the outskirts (central panels) and in groups (right panel for EDisCS) and in the field (right 
panel for ICBS). Black vertical dashed lines represent the range in stellar mass probed by the surveys. Curves are not normalized. 
The blue ICBS cluster mass function is omitted due to poor number statistics. 



- Galaxies in cluster (R/R200 ^ 1), in groups, and in the field 
seem to follow similar mass distributions, because no statis- 
tical differences can be detected. Hence, we find no evidence 
for a dependence of the shape of the galaxy stellar mass func- 
tion on global environment at z = 0.3 - 0.8. 

- By comparing of the ICBS mass functions with mass func- 
tions in the local Universe ( Vulcani et al.ll20ll] for clusters 
and Calvi et al. in preparation for the field), we have found 
that the evolution of the shape of the mass function from 
z ~ 0.4 to z ~ 0.07 is the same in the field and in clusters, 
hence it turns out to be independent on global environment. 
However, we stress that even though the shape of the stellar 
mass function from high to low redshift in clusters and in the 
field is similar, the amount of growth as a function of cosmic 
time is not the same at the low- and high-mass end. 

- Virialized regions of clusters at various clustercentric dis- 
tances present very similar mass functions, as it happens also 
for galaxies within and outside R/R200 - 1. 

- Subdividing galaxies in terms of color, our results suggest 
that in clusters, groups and field, red and blue galaxies are 
regulated by different mass functions. When comparing the 
mass function in different environments separately for blue 
and red galaxies, no differences are detected. 

All ou r results have bee n confirmed also by the analysis of 
the best-fit [Schechterl d 19761) parameters. 

To summarize, our results show that global environment does 
not seem to have a prominent role in shaping galaxy stella r mass 
functions. In contrast, as presented inl Vulcani et al.l(l2012l) . local 
environment seems to have a different influence in determining 
the mass distribution. This suggests that the most fundamental 
of all galaxy properties, the galaxy mass, is not much dependent 



of cluster mass, but do depend on local scale processes. Global 
and local environment are clearly related to different physical 
processes, and their different role in altering galaxy properties is 
important to understanding the drivers of galaxy formation and 
evolution. 



Appendix A: Is the mass function simply a mirror of 
the luminosity function? 

Using exactly the same samples of 95.11 in this Appendix we 
build the luminosity functions, to test whether luminosity and 
mass functions give us the same results and therefore the same 
information. Both for ICBS and EDisCS, we take into account 
the absolute magnitude My, derived as presented in 93.11 and 
EH 

In the same way we built mass functions, in each magni- 
tude bin, we sum all galaxies belonging to the environment under 
consideration to obtain the total number of galaxies. The width 
of each magnitude bin is 0.4 dex. In each plot, histograms are 
normalized using the total integrated luminosity. 

In ICBS (left panel of Figure IATTT i. we compare cluster re- 
gions and the field. The two distributions are different: in clusters 
the number of more luminous galaxies is proportionally higher 
than in the field. The K-S test supports our finding: giving a prob- 
ability of ~ 1.4% it excludes the similarity of the distributions. 

In EDisCS (right panel of Figure lA.lb . we compare cluster 
regions and groups. Again, both from a visual inspection of the 
plot and from the K-S test (Pk-s ~ 2.5%), we can conclude that 
the two distributions are different. Groups seem to have a higher 
number of more luminous galaxies than in clusters. 
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Fig. A.l. Luminosity function in the different environments. Left panel: ICBS cluster regions (black crosses) and field (red filled 
squares). Right panel: EDisCS clusters (black crosses) and groups (red filled squares). In each panel, luminosity functions are 
normalized using the total integrated luminosity. At the bottom of each panel, the K-S probabifities are given. Significant differences 
are evident between the luminosity functions in different environments. 



Therefore, both using ICBS and EDisCS data, the luminos- 
ity function of galaxies in different global environments is statis- 
tically different, while the mass functions are indistinguishable 
(see 35.11 ). This shows that studying the luminosity function does 
not give direct information about the mass function. This is be- 
cause galaxies do not have all the same colors, hence one single 
mass-to-light ratio, such as that of passively evolving galaxies. 
As a consequence, when we derive masses from luminosities we 
are not simply multiplying luminosities by a constant factor 

Acknowledgements. We thank the anonymous referee for her/his detailed and 
careful comments and suggestions which helped us to improve the manuscript. 
We thank Lucia Pozzetti for providing us z-COSMOS data for the field mass 
functions and useful discussion. We also thank Micol Bolzonella for for her 
suggestions. BV and BMP acknowledge financial support from ASl contract 
1/016/07/0 and ASl-lNAF 1/009/10/0. BV also acknowledges financial support 
from the Fondazione Ing. Aldo Gini. GDL acknowledges financial support 
from the European Research Council under the European Community's Seventh 
Framework Programme (FP7/2007-2013)/ERC grant agreement n. 202781. 



References 

Adami C, Biviano A., Mazure A., 1998, A&A, 331, 439 

Aragon-Salamanca, A. et al. in preparation 

Baldry, I. K., et al. 2004, ApJ, 600, 681 

Baldly, I. K., et al. 2006, MNRAS, 373, 469 

Baldry, I. K., et al. 2008, MNRAS, 388, 945 

Balogh M., Navarro J., Morris S., 2000, gfc.conf 

Baldly, I. K., Driver S. P, Loveday, J., et al. 2012, MNRAS, 421, 621 

Balogh M. L., Chiistlein D., Zabludoft' A. I., Zaiitsky D., 2001, ApJ, 557, 1 17 

Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 

Bekki, K. 2009, MNRAS, 399, 2221 

Bell, E. R, & de Jong, R. S. 2001, ApJ, 550, 212 

Bell, E. F., Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289 

Bell, E. R, et al. 2004, ApJ, 608, 752 

Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 

Bolzonella M., Miralles J.-M., Pello R., 2000, A&A, 363, 476 

Bolzonella M., et al., 2010, A&A, 524, A76 

Borch, A., et al. 2006, A&A, 453, 869 

Boselli, A., & Gavazzi, G. 2006, PASP 118, 517 

Brammer G. B., Whitaker, K. E., van Dokkum, P G., et al. 201 1, ApJ, 739, 24 
Brunner R. J., & Lubin, L. M. 2000, AJ, 120, 2851 



Bundy, K.,et al. 2005, ApJ, 625, 621 
Bundy, K., et al. 2006, ApJ, 651, 120 
Byrd, G., & Valtonen, M. 1990, ApJ, 350, 89 

Calvi R., Poggianti B. M., Vulcani B., 2011, arXiv. larXiv: 1105.3683 1 
Cappellaii, M., McDermid, R. M., Alatalo, K., et al. 2012, Nature, 484, 485 
Caputi, K. I., Cirasuolo, M., Dunlop, J. S., et al. 2011, MNRAS, 413, 162 
Conroy, C, White, M., & Gunn, J. E. 2010, ApJ, 708, 58 
De Lucia, G., et al. 2004, ApJ, 610, L77 

De Lucia G., et al., 2007, MNRAS, 374, 809 

Dekel G., Farberovich V, Fleurov V, Sofl'er A., 2009, arXiv. larXiv:0911. 15371 

Desai, V, et al. 2007, ApJ, 660, 1 151 

Dressier, A., et al. 2004, ApJ, 617, 867 

Driver S. P, et al., 2006, MNRAS, 368, 414 

Drory, N., et al. 2005, ApJ, 619, L131 

Drory N., Fisher D. B., 2007, ApJ, 664, 640 

Drory, N., Bundy, K., Leauthaud, A., et al. 2009, ApJ, 707, 1595 

Dutton, A. A., Conroy, C, van den Bosch, R C, et al. 201 1, MNRAS, 416, 322 

Eisner F, Feulner G., & Hopp, U. 2008, A&A, 477, 503 

Fasano, G., et al. 2006, A&A, 445, 805 

Fasano G., et al., 2010, yCat, 740, 41490 

Fontana, A., et al. 2004, A&A, 424, 23 

Fontana, A., et al. 2006, A&A, 459, 745 

Gehrels, N. 1986, ApJ, 303, 336 

Gladders M. D., Yee H. K. C, 2000, AJ, 120, 2148 

Gonzalez, A. H., et al. 2001, ApJS, 137, 117 

Gonzalez, V, Labbe, 1., Bouwens, R. J., et al. 2011, ApJ, 735, L34 

Guo Y., et al., 2009, MNRAS, 398, 1 129 

Gunawardhana M. L. P, et al., 201 1, MNRAS, 415, 1647 

Gunn, J. E., & Gott, J. R., Ill 1972, ApJ, 176, 1 

Gwyn, S. D. J., & Hartwick, R D. A. 2005, AJ, 130, 1337 

HalHday, C, et al. 2004, A&A, 427, 397 

Holden B. R, et al., 2009, ApJ, 693, 617 

Kajisawa, M., Ichikawa, T., Tanaka, I., et al. 2009, ApJ, 702, 1393 

Kauftmann G., et al., 2003, MNRAS, 341, 54 

Kovac K., et al., 2010, ApJ, 718, 86 

Kroupa, R 2001, MNRAS, 322, 231 

llbert, O., et al. 2010, ApJ, 709, 644 

lovino A., et al., 2010, A&A, 509, A40 

Larson R. B., Tinsley B. M., Caldwell C. N., 1980, ApJ, 237, 692 

Marchesini D., van Dokkum P. G., Forster Schreiber N. M., Franx M., Labbe 1., 

Wuyts S., 2009, ApJ, 701, 1765 
Mercurio A., Haines C. P., Merluzzi P., Busaiello G., Smith R. J., Raychaudhuiy 

S., Smith G. P, 2010, jena.conf 
Mihos J. C, Hemquist L., 1996, ApJ, 464, 641 



18 



Benedetta Vulcani et al.: No dependance of the mass function on global environment 



Milvang-Jensen, B., et al. 2008, A&A, 482, 419 
Mo H. J., White S. D. M., 1996, MNRAS, 282, 347 
Moore, B., et al. 1996, Nature, 379, 613 

Mortlock, A., Conselice, C. J., Bluck, A. F. L., et al. 2011, MNRAS, 413, 2845 
Moster, B. R, Somerville, R. S., Maulbetsch, C, van den Bosch, F. C, Maccio, 

A. v., Naab, T., & Oser, L. 2010, ApJ, 710, 903 

Muldrew, S. I., Croton, D. J., Skibba, R. A., et al. 201 l. larXiv:1109T6328] 
Muzzin, A., Marchesini, D., van Dokkum, R G., et al. 2009, ApJ, 701, 1839 
Pasquali A., et al., 2008, ApJ, 687, 1004 

Pasquali, A., van den Bosch, F. C, Mo, H. J., Yang, X., & Somerville, R. 2009, 

MNRAS, 394, 38 
Peng, Y.-j., et al. 2010, ApJ, 721, 193 
Pello, R., et al. 2009, A&A, 508, 1 173 
Poggianti, B. M., et al. 2006, ApJ, 642, 188 
Poggianti, B. M., et al. 2008, ApJ, 684, 888 
Pozzetti, L., et al. 2007, A&A, 474, 443 
Pozzetti L., et al., 2010, A&A, 523, A13 
Rudnick, G., et al. 2001, AJ, 122, 2205 
Rudnick, G., et al. 2003, ApJ, 599, 847 
Rudnick, G., et al. 2006, ApJ, 650, 624 
Rudnick, G., et al. 2009, ApJ, 700, 1559 
Salpeter, E. E. 1955, ApJ, 121, 161 

Sandage, A., Tammann, G. A., & Yahil, A. 1979, ApJ, 232, 352 

Schechter R, 1976, ApJ, 203, 297 

Sheth R. K., Tormen G., 2002, MNRAS, 329, 61 

Schawinski K., Virani S., Simmons B., Urry C. M., Treister E., Kaviraj S., 

Kushkuley B., 2009, ApJ, 692, L19 
Spmiello, C, Trager, S. C, Koopmans, L. V. E., & Chen, Y. R 2012, ApJ, 753, 

L32 

Taylor, E. N., et al. 2009, ApJS, 183, 295 

Thomas D., Maraston C, Schawinski K., Sarzi M., Silk J., 2010, MNRAS, 404, 
1775 

Thomas, J., Saglia, R. R, Bender, R., et al. 201 1, MNRAS, 415, 545 

Treu, T, Gavazzi, R., Gorecki, A., et al. 2009, ApJ, 690, 670 

van den Bosch F. C., Aquino D., Yang X., Mo H. J., Pasquah A., Mcintosh D. H., 

Weinmann S. M., Kang X., 2008, MNRAS, 387, 79 
van den Bosch F. C., Pasquali A., Yang X., Mo H. J., Weinmann S., Mcintosh 

D. H., Aquino D., 2008, arXiv, arXiv:0805.0002 
von der Linden A., Wild V., Kauffmann G., White S. D. M., Weinmann S., 2010, 

MNRAS, 404, 1231 
van Dokkum, P G., & Conroy, C. 201 1, ApJ, 735, L13 
Vulcani, B., et al. 2010, ApJ, 710, LI 
Vulcani, B., et al. 2011, MNRAS, 412, 246 

Vulcani, B., Poggianti, B. M., Fasano, G., et al. 2012, MNRAS, 420, 1481 
Weinmann S. M., van den Bosch F. C., Yang X., Mo H. J., 2006, MNRAS, 366, 
2 

White, S. D. M., et al. 2005, A&A, 444, 365 

Wild v., Kauffmann G., Heckman T., Chariot S., Lemson G., Brinchmann J., 
Reichard T., Pasquali A., 2007, MNRAS, 381, 543 



19 



